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We study theoretically the shift of the superconducting transition temperature (T c ) 
under uniaxial compression in /3-type organic superconductors, /3-(BEDT-TTF) 2 l3 and 
/3-(BDA-TTP)2X[X=SbF6,AsFg], in order to clarify the electron correlation, the spin 
frustration, and the effect of dimerization. The transfer integrals are calculated by 
the extended Hiickel method assuming the uniaxial strain, and the superconducting 
state mediated by the spin fluctuation is solved using Eliashberg's equation with the 
fluctuation-exchange approximation. The calculation is carried out on both the dimer- 
ized (one-band) and nondimerized (two-band) Hubbard models. We have found that (i) 
the behavior of T c in /3-(BEDT-TTF)2i3 with a stronger dimerization is well reproduced 
by the dimer model, while that in weakly dimerized /3-BDA-TTP salts is rather well re- 
produced by the two-band model, and (ii) the competition between the spin frustration 
and the effect induced by the fluctuation is important in these materials, which causes 
the nonmonotonic shift of T c against uniaxial compression. 

KEYWORDS: superconductivity, organic superconductor, uniaxial compression, spin frustration, spin fluc- 
tuation 

1. Introduction 

The superconductivity of the K-type BEDT-TTF salts has fascinated experimental 
and theoretical researchers with a high transition temperature T c over 10 K. 1 ) In the 
K-type BEDT-TTF salts, donors are considered to form dimers. 2 ^ Then, the split anti- 
bonding HOMO band is half-filled with the effective on-site Coulomb interaction (U) 
close to the bandwidth. It is considered to be a strongly correlated electron system, 
where the superconductivity is mediated by the antiferromagnetic spin fluctuation. 3 ~ 7 ) 
The spin fluctuation is suppressed by the geometrical spin frustration since the donors 
locate on an anisotropic triangular lattice. The effect becomes strongest on a regular 
triangular lattice, which is considered to be realized in k-(BEDT-TTF)2Cu2(CN)3. 8,9 ) 
The organic superconductors with /3-type donor arrangement also exhibit a relatively 
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high T c reaching 8 K. As shown in Fig. 1, donors stack unidirectionally and dimerize 
on a triangular lattice. However, the dimerization of the molecules is weaker than that 
of the K-type salts. Then, as an effective model for the /3-type salt, not only the dimer 
model but also the nondimerized original two-band model should be taken into account. 

Uniaxial compression is a powerful technique for studying the strongly correlated 
electron systems since transfer integrals can be controlled selectively 10 -* The spin frus- 
tration and the electron correlation are controlled by uniaxial compression. In fact, 
the shift of T c under the uniaxial compression in k-(BEDT-TTF) 2 Cu(NCS) 2 can be 
explained by the term of the spin frustration on an anisotropic triangular lattice com- 
posed of dimerized molecules. 11 ) Our group reported the shift of T c in /3-type salts, 
/3-(BDA-TTP) 2 X [X=SbF 6 , AsF 6 ], under uniaxial compression parallel and perpendic- 
ular to the molecular stack direction. 12 ) Under the compression perpendicular to the 
stack, T c shows a nonmonotonic behavior, once increases and reaches a maximum and 
decreases subsequently, as a function of the applied piston pressure. On the basis of 
the dimer model of the anisotropic triangular lattice, this behavior has been attributed 
to the result of the competition of the enhancements of the spin fluctuation and the 
spin frustration effect. However, the data for /3- (BDA-TTP ) 2 AsF 6 was not well repro- 
duced. Furthermore, for /3-(BEDT-TTF) 2 I 3 , T c exhibits nonmonotonic behaviors for 
both compression directions, in strong contrast to that of the BDA-TTP salts. 13 ^ 

In order to clarify the behavior of T c under uniaxial compression and obtain insight 
into the superconductivity of /3-type salts, we study theoretically the superconductivity 
mediated by the spin fluctuation on both the original two-band Hubbard model and the 
dimerized Hubbard model using the fluctuation exchange (FLEX) approximation 14,15 ) 
for /3-(BEDT-TTF) 2 I 3 and 0- (BDA-TTP) 2 X[X=SbF 6 ,AsF 6 ]. Measurements of STM 
and specific heat are consistent with the rf-wave symmetry of BDA-TTP salts. 16 ' 17 * 1 It 
is noted that the zero bias conductance peak due to the Andreev bound state 18 ' 19 ^ has 
been reported very recently. 20 ) Thus, it is natural to employ the spin fluctuation theory, 
which can realize a <i-wave gap function beyond the BCS theory. 

Although some of the results have been reported in our previous papers, 12 ' 13,21 - ) in 
this paper, we note the difference in the strength of the dimerization among the three 
materials, which is evaluated by the ratio of the transfer integrals along the stacking 
direction, namely, t p \ and t p2 , as shown in Fig. 1. We have found that the dimer model 
is suitable to the strongly dimerized /3-(BEDT-TTF) 2 l3 among the three materials, 
tpi/t P 2 = 2.99, while the original two-band model is suitable for the weakly dimerized 
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0- (BDA-TTP)2AsFg among the three materials, t p \/t P 2 = 2.02. We have also found that 
the nonmonotonic behavior of T c is due to the competition between the spin frustration 
and the effect of the self-energy induced by the spin fluctuation. 




Fig. 1. (Color online) Lattice structure on a two-band model (left) and the dimer model (right) for 
/3-type salts. Donor molecules locate at A and B sites, and form dimers, as depicted by the dotted 
ellipsoids. 



2. Formulation 

In the following, we introduce the models that we employ in this study. Putting a 
donor molecule into a site on a lattice, the crystal structure of /3-type salts is modeled 
by a triangular lattice with A and B sublattices, as shown in Fig. 1. Here, we distort the 
lattice in order to make the reciprocal lattice rectangular. Then, the direction perpen- 
dicular to the donor stack direction is distorted, as illustrated in Fig. 1. However, this 
distortion makes no changes in physical quantities, such as T c , because the topology of 
the lattice is invariant. According to the dimer model, the sublattices are dimerized, 
and the original two-band model is reduced to the one-band model (dimer model), as 
shown in Fig. 1. In the dimer model, hopping integrals are given by t x = —t^jl + t c , 
t xy = —t q i/2, and t y = —t P 2/2 using the hopping integrals in the original two-band 
model. 22 ** Hopping integrals under uniaxial compression are also estimated by the ex- 
tended Hiickel method 23 ^ assuming uniform displacements of molecules from the struc- 
ture at ambient pressure. 12 ^ The validity of the Huckel method as a first approximation 
is implied because it can predict the round Fermi surface, which is consistent with 
the Shubnikov-de Haas (SdH) or de Haas- van Alphen (dHvA) measurements. 24 ' 25 ^ The 
uniaxial strain under uniaxial compression is assumed to be 0.3%/kbar according to 
the results of X-ray measurements. 26 ^ We also assume that the encapsulation of the 
sample with epoxy, which is prepared for uniaxial compression, induces an isotropic 
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Table I. Transfer integrals under uniaxial compression in the /3-type salts, where unit of transfer 
integrals is 10 _2 eV. 





Pressure [kbar] 


to 


t p i 


tp2 


tql 


tq2 


/3-(BEDT-TTF) 2 I 3 





-5.22 


24.5 


8.23 


12.6 


6.91 


//stack 5 


-4.75 


26.8 


9.16 


13.0 


7.07 


_L stack 5 


-5.98 


24.0 


8.15 


13.1 


7.25 


/3-(BDA-TTP) 2 SbF 6 





-0.45 


15.0 


6.43 


8.34 


9.04 


//stack 5 


-0.02 


16.3 


7.12 


8.58 


9.23 


_L stack 5 


-0.85 


14.8 


6.38 


8.68 


9.60 


/3-(BDA-TTP) 2 AsF 6 





-0.36 


15.5 


7.74 


8.52 


9.79 


/ / stack 5 


0.07 


16.8 


8.45 


8.72 


9.99 


_L stack 5 


-0.78 


15.2 


7.69 


8.83 


10.5 



pressure of lkbar in the sample. We show the obtained transfer integrals under uniaxial 
compression in /3-type salts in Table I. 

In this paper, we employ both the two-band and one-band Hubbard models with 
the repulsion U on a triangular lattice, 
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in standard notations and calculate T c by solving Eliashberg's equation within the FLEX 
approximation. 27 ' 28 ^ In the two-band model, Green's function G, the self-energy S, the 
susceptibility \i an d the gap function are all 2 x 2 matrices, such as G a p{k), where 
a,/3=A, B and k = (k,iu n ) with u n = (2n — 1)ttT being the Matsubara frequency for 
fermions. 

The self-energy is given by 



S^(*) = ^X)G ! a/J (fc- 9)^(9), 
where the FLEX effective interaction V(q) for the self-energy is 
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with the irreducible susceptibility 

7&(q) = ~ E G ^ k + l) G Uk)- (4) 

Here, we denote q = (q, ie{) with e\ = 2nlT being the Matsubara frequency for bosons, 
and N as the number of fc-points on a mesh. 
With Dyson's equation, 

Hk)~ 1 ] afi = [G^)- 1 ]^-^^), (5) 

where G° is the bare Green's function G^Jk) = [(iu n + fi — with e® k being the 

bare energy, we have solved eqs. (2)- (5) self-consistently to obtain Green's functions 
under FLEX interactions. 

Then, T c is obtained from Eliashberg's equation for the spin-singlet pairing, 

rp A,B 

k' a',/3' 

y.G aa ,(k , )Gw(-k , )<j> Q , fi ,(k'), (6) 

where is the gap function, and the spin singlet pairing interaction V sm (q) is given as 

+U5 af) . 

From eq. (6), T c is determined as the temperature at which the maximum eigenvalue 
A becomes unity. However, because the performance of the computer is not sufficient 
for deducing T c directly, we evaluate A at a constant temperature of T/t q i = 0.01 as 
a function of the piston pressure for the uniaxial compression. The larger value of A 
corresponds to the higher T c . 

At this temperature, the fc-dependence of the spin susceptibility x s is expressed by 
the diagonalized component 

2 

+ IXab! 2 , (8) 
where the matrix of the spin susceptibility is 

(9) 

a/3 

Throughout this study, we consider t q i under ambient pressure as a unit of energy, U = 
10 on the two-band model and U = 2t p \ on the dimer model. The system with U — 10 
corresponds to the strongly correlated electron system. In the numerical calculation, we 
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Fig. 2. (Color online) Experimental data of T c under uniaxial compression in /3-(BEDT-TTF)2l3 
defined by the midpoint of the transition. 13 ) 

consider N = 128 x 128 fc-point meshes and the Matsubara frequencies from e n from 
-(2N C - 1)ttT to {2N C - 1)ttT with N c = 4096. 

3. Result 

3.1 I3-(BEDT-TTF) 2 I 3 

We begin with /3-(BEDT-TTF) 2 I 3 in which the dimerization is the strongest among 
the three salts. Although it is known that the T c of /3-(BEDT-TTF) 2 l3 jumps from 1.5K 
at ambient pressure up to 8K at a pressure of lkbar, 29 ) only the 8K superconductivity 
phase is observed under epoxy encapsulation because of the slight pressure induced by 
the epoxy. 13 ** Figure 2 shows that the experimental T c increases with uniaxial compres- 
sion along both the / / stack and _Lstack directions up to 4kbar, and decreases in further 
piston pressures. 13 - ) It is very interesting that T c shows a nonmonotonic behavior against 
the uniaxial compression. 




Fig. 3. (a) Fermi surface (solid line) and nodes of gap function (dotted line) and (b) \s obtained 
using the two-band model for /3-(BEDT-TTF)2l3 at ambient pressure. 

We first examine the gap function and the Fermi surface. The gap function obtained 
using eq. (6) is a d-wave that changes sign four times on the Fermi surface, as shown 
in Fig. 3. The Fermi surface obtained using the two-band model is closed owing to the 
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periodicity of the Brillouin zone (BZ), which is consistent with the previous result 30 ) 
and dHvA measurements. 24, 25 - ) The Fermi surface in the two-band model is almost the 
same as that in the dimer model. The peak of Xs is observed at the wave number (jr, 0). 
The gap function changes sign along the nesting vector Q = (vr,0). We consider that 
/3-(BEDT-TTF) 2 l3 is close to the SDW phase because the peak value of Xs is large 
(~ 6). 




Fig. 4. (Color online) A under uniaxial compression on the dimer model (left panel) and two-band 
model (right panel) for /3-(BEDT-TTF) 2 I 3 . 



Next, we see the shift of A, i.e., T c under uniaxial compression. Assuming a linear 
response to the uniaxial strain, the experimental results show an increase in T c under low 
piston pressures along both directions. Here, the absolute value of A is not important. 
Thus, we focus on the shift of A around ambient pressure. In the case of the dimer model, 
A increases with uniaxial compression along both the / / stack and _Lstack directions at 
low pressures, as shown in Fig. 4. The trend of the increasing A along both compression 
directions agrees with the experimental result at low pressures, although the larger A 
under compression perpendicular to the stacks does not match with the order of the 
experimental T c . 

On the other hand, in the case of the two-band model, A increases with uniaxial 
compression along the / / stack direction, while A decreases along the _Lstack direction. 
This result does not match the experimental results. Thus, the dimer model is more 
suitable for /3-(BEDT-TTF)2l3, which reproduces the shift of T c at low pressures, than 
the two-band model. However, the nonmonotonic behavior of T c , i.e., the subsequent 
decrease in T c at high pressures, can be reproduced by neither model. 
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Pressure [kbar] 

Fig. 5. (Color online) Experimental data of T c under uniaxial compression in /3-(BDA-TTP)2SbF6 
denned by the midpoint of the transition. 12 ) 

3.2 f3-(BDA-TTP) 2 SbF 6 

Next, we move on to /3-(BDA-TTP)2SbFg. Under uniaxial compression along the 
/ /stack direction, experimental data of T c increases with pressure, reaching a maximum 
at 3kbar, and decreases in further pressures, as shown in Fig. 5. On the other hand, 
along the _Lstack direction, T c decreases monotonically with pressure. 




Fig. 6. (a) Fermi surface (solid line) and nodes of gap function (dotted line) and (b) \s obtained 
using the two-band model for /3-(BDA-TTP) 2 SbF 6 at ambient pressure. 

The gap function is also calculated to be a cf-wave, as shown in Fig. 6(a). The 
obtained Fermi surface and symmetry of the gap function are consistent with previous 
calculations 31 ' 32 ^ and STM measurements, 16 - ) indicating the four-fold symmetry of the 
gap function. The (i-wave is mediated by the spin fluctuation induced by the nesting 
vector. Thus, it is also reasonable to apply the spin fluctuation theory to this system. 
As shown in Fig. 6(b), the peak of Xs is found at the wave number (^j7r,|7r); however, 
the peak value of Xs is small (~ 1.4). 

Then, we observe the shift of A under uniaxial compression in Fig. 7. In the case of 
the dimer model, A increases with pressure reaching a maximum at 2kbar, and decreases 
in further pressures under uniaxial compression along the / / stack direction, while along 
the _Lstack direction, A decreases monotonically with increasing pressure. This behavior 
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dimer model two-band model 




Pressure[kbar] Pressure[kbar] 

Fig. 7. (Color online) A under uniaxial compression on the dimer model (left panel) and the two- 
band model (right panel) for /3- (BDA-TTP) 2 SbF 6 . 

is consistent with the experimental results. On the other hand, in the case of the two- 
band model, A increases monotonically with pressure under uniaxial compression along 
the //stack direction, while A decreases monotonically under uniaxial compression along 
the _Lstack direction. 

In /3-(BDA-TTP) 2 SbF 6 , a maximum A value is obtained on the dimer model, which 
is not obtained in /3-(BEDT-TTF) 2 l3. In the following, we discuss the origin of the 
nonmonotonic behavior of A under uniaxial compression calculated in the dimer model. 

First, we focus on {U 2 xs)b that is the mean value of U 2 Xs within the BZ, which 
is a significant term in Eliashberg's equation (6). In the case of the dimer model, the 
on-site Coulomb interaction U = 2t pl depends on the uniaxial compression. Thus, the 
value of U 2 in the effective interaction is also important when we discuss the uniaxial 
compression dependence of A. As shown in Fig. 8, under uniaxial compression along 




Pressure [kbar] 

Fig. 8. (Color online) {U 2 xs)b under uniaxial compression on the dimer model for /3-(BDA- 
TTP) 2 SbF 6 . 

the //stack direction, (U 2 x s )b increases monotonically with pressure, while {U 2 x s )b 
decreases under uniaxial compression along the _Lstack direction. Here, we note the 
geometrical spin frustration induced by the triangular lattice. We denote t'(= t y ) as the 
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hopping integral along the dimerized direction and t(= (t x + t xy )/2) as the mean value 
of other hopping integrals in the dimer model. In this material, t'/t < 1 is satisfied 
at ambient pressure, and then the geometrical spin frustration becomes strong under 
uniaxial compression along the / / stack direction, while it becomes weak under uniaxial 
compression along the _Lstack direction, which is consistent with the results of (U 2 Xs)b- 

3.8 
3.7 
3.6 

CO 

A 3 - 5 

<N_ 3.4 
^3.3 
V 3.2 
3.1 
3 
2.9 

Fig. 9. (Color online) (|G| 2 )b under uniaxial compression on the dimer model for /3-(BDA- 
TTP) 2 SbF 6 . 

Next, we show in Fig. 9 the mean value of \G\ 2 in the BZ, which is another im- 
portant term in eq. (6). Under uniaxial compression along the //stack direction, the 
value of (|G| 2 )b decreases monotonically with increasing pressure, while that increases 
monotonically along the _Lstack direction. Under uniaxial compression along the / / stack 
direction, the increase in (U 2 x s )b strongly competes with the decrease in (|G| 2 ) B . Thus, 
the nonmonotonic behavior reaching a maximum at the value of A is considered to be 
caused by such competition. 

0.22 
0.2 
030.18 

A 

W 0.16 
m 0.14 

V 

0.12 
0.1 
0.08 

1 2 3 4 5 6 

Pressurefkbar] 

Fig. 10. (Color online) (|im£|)B under uniaxial compression on the dimer model for /3-(BDA- 
TTP) 2 SbF 6 . 

Now, we consider the origin of the uniaxial compression dependence in (|C| 2 )b- G 
is suppressed by the imaginary part of self-energy ImE. In Fig. 10, we show the com- 
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pression dependence of (|ImS|)B, which is the mean value of |ImS| in the BZ. Under 
uniaxial compression along the / /stack direction, the value of (|ImS|)B increases mono- 
tonically with pressure, while it decreases monotonically along the _Lstack direction. 
These results are consistent with the results of (|G| 2 )b- Thus, under uniaxial compres- 
sion along the / /stack direction, the increase in the value of (|Ihi£|)b reduces the value 
of (|G| 2 )b, which competes with the increase in the value of (U 2 x s )b 

Compared with the experimental results, the low-pressure behavior of T c can be 
reproduced by both the dimer and two-band models. We observe that the dimer model is 
more suitable than the two-band model, because the maximum of T c can be reproduced 
by the dimer model. 

3.3 (3-(BDA-TTP) 2 AsF 6 

8 
7.5 

7 

g 6.5 
K° 6 

5.5 
5 

4.5 

1 2 3 4 5 6 
Pressure [kbar] 

Fig. 11. (Color online) Experimental data of T c under uniaxial compression in /3-(BDA-TTP)2AsF6 
defined by the midpoint of the transition. 12 ' 

Finally, we move on to /3-(BDA-TTP) 2 AsF 6 . As shown in Fig. 11, under uniaxial 
compression along the / / stack direction, the experimental data of T c is almost con- 
stant up to 3kbar and increases with pressure, reaching a maximum at 4.5kbar, and 
decreases in further pressures. On the other hand, along _Lstack direction, T c decreases 
monotonically with increasing pressure. 12 - ) 

As shown in Fig. 12, the calculated gap function is also a d-wave, and the shape 
of the Fermi surface is consistent with previous results 31 ^ and STM measurements. 16 -* 
The peak of Xs is found at the wave number (7r,|7r), but the peak value of Xs is small 
(~ 1.2). We can assume that this material is far from the SDW phase owing to the 
small value of Xs- 

Figure 13 shows that in the dimer model, A increases only very slightly up to 1 
kbar and decreases monotonically with increasing pressure under uniaxial compression 
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Fig. 12. (a) Fermi surface (solid line) and nodes of gap function (dotted line) and (b) %s obtained 
using the two-band model in /3-(BDA-TTP) 2 AsF 6 at ambient pressure. 



dimer model 



two-band model 




2 3 4 5 

Pressure [kbar] 



2 3 4 5 

Pressure[kbar] 



Fig. 13. (Color online) A under uniaxial compression on the dimer model (left panel) and two-band 
model (right panel) for /3-(BDA-TTP) 2 AsF 6 . 



along the / / stack direction and A also decreases monotonically along _Lstack direction, 
which does not fit with the experimental results. 12 - ) On the other hand, in the two- 
band model, A increases with pressure, reaching a maximum at 3.5kbar, and decreases 
in further pressures under uniaxial compression along the //stack direction, while A 
decreases monotonically along the _Lstack direction. These behaviors of T c in the two- 
band model are consistent with those observed by experiment. Thus, the two-band 
model is more suitable for reproducing the behavior of T c in /3-(BDA-TTP) 2 AsF 6 . 




39 1 1 1 1 1 1 1 
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Fig. 14. (Color online) (U 2 xs)b under uniaxial compression on the two-band model for /3-(BDA- 
TTP) 2 AsF 6 . 
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Fig. 15. (Color online) (|G| 2 )b under uniaxial compression on the two-band model for /3-(BDA- 
TTP) 2 AsF 6 . 

Hereafter, we discuss the origin of the nonmonotonic behavior of A in the two-band 
model. As shown in Fig. 14, in the two-band model, {U 2 xs)b decreases with increasing 
pressure under uniaxial compression along both the / / stack and _Lstack directions. In 
Fig. 15, we show the result of (|C| 2 )b that is the mean value of the largest eigenvalue of 
\G\ 2 in the BZ. Under uniaxial compression along both the //stack and _Lstack direc- 
tions, the value of (|C| 2 )b increases monotonically with pressure. Thus, under uniaxial 
compression along the //stack direction, the decrease in (U 2 x s )b and the increase in 
(|G| 2 )b strongly compete, as in the case of /3-(BDA-TTP) 2 SbF 6 . The strong competi- 
tion between (U 2 x s )b and (|G| 2 )b is the origin of the nonmonotonic behavior of the 
value of A. 

0.12 
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CQ 

s 

^ 0.105 
D.I 
0.095 
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Fig. 16. (ImE)B under uniaxial compression on the two-band model for /3-(BDA-TTP)2AsFg. 

The uniaxial compression dependence of (|G| 2 )b is explained by the largest eigen- 
value of (|Ihi£|)b as shown in Fig. 16, as in the case of /3-(BDA-TTP) 2 SbF 6 . Under uni- 
axial compression along the both //stack and _Lstack directions, the value of (|ImS|)B 
decreases monotonically with increasing pressure. In particular, along the / / stack di- 
rection, the value of (|Im£|)B decreases rapidly. This decrease in the value of (|Im£|)B 
increases the value of (|G| 2 )b. 
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Table II. Relationship between the dimerization factor (t p i/f P 2) and the suitable model for the 
three /3-type salts. 





^-(BEDT-TTF) 2 I 3 


^-(BDA-TTP) 2 SbF 6 


/?-(BDA-TTP) 2 AsF 6 


dimerization factor (t p i/i P 2) 


2.98 


2.30 


2.00 


dimer model 








X 


2 band model 


X 









4. Discussion 

As a summary of the results, we show Table II that indicates the suitable model 
for explaining the experimental results of T c , mainly focusing on the reproducibility of 
the experimental T c in a low pressure range, where the linearity of the uniform strain 
seems to be a good approximation. The strength of dimerization is estimated by the 
dimerization factor (t p i/t p2 ) because the direction of dimerization corresponds to the 
direction of t p i (intradimer) and t P 2 (interdimer). 

Table II shows the following results, (i) In the case of /3-(BEDT-TTF) 2 I 3 , the low- 
pressure behavior of T c cannot be reproduced by the two-band model, but it can be 
reproduced by the dimer model, (ii) In the case of /3-(BDA-TTP) 2 SbF 6 , the low-pressure 
behavior of T c can be reproduced by both the dimer model and the two-band model. We 
also observe that the dimer model is more suitable than the two-band model, because 
the maximum of T c can be reproduced by the dimer model, (iii) In the case of /3-(BDA- 
TTP)2AsFg, the low-pressure behavior of T c cannot be reproduced by the dimer model, 
but it can be reproduced by the two-band model. We also observe that the two-band 
model can reproduce the maximum of T c . Thus, when the dimerization factor (t p i/t P 2) 
is large, which corresponds to the case of strong dimerization, experimental results are 
well-reproduced by the dimer model. On the other hand, when the dimerization factor 
is small, the results of the two-band model are suitable for reproducing the experimental 
results. 

One may consider that the results of the dimer model are reproduced by the two- 
band model. However, we find that the results of the dimer model are not always the 
same as those of the two-band model with finite U = I0t ql because the dimer model 
corresponds to the two-band model with U — > oo. The dimer model has the effective 
t^dimer = 2t p i at the limit of U — > oo. In fact, we cannot find the value of U, using which 
the results of /3-(BEDT-TTF) 2 i3 in the dimer model are reproduced in the two-band 
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model. 

Our calculation cannot reproduce the experimental results at high pressures. Our 
assumption of the uniform displacements of molecules at high pressures may be in- 
correct. Taking account of the nonlinear distortion of the lattice against the uniaxial 
compression is one of the future problems. Thus, in this paper, our calculation should 
be compared with the experimental results obtained at low pressures. 

As for the difference among the three materials, we observe that the calculated 
Fermi surfaces of /3-(BDA-TTP) 2 SbF 6 and /3-(BDA-TTP) 2 AsF 6 are slightly thinner 
than that of /3-(BEDT-TTF) 2 I 3 , which is consistent with the recent experimental results 
of magnetoresistance. 33 ' 34 ) 

The structure of the spin susceptibility, which is very sensitive to the lattice struc- 
ture, is important for the spin-fluctuation-mediated superconductivity. Thus, it may be 
essential that we employ the hopping integrals obtained by the first-principles calcula- 
tion in order to discuss the effect of the uniaxial compression in more detail. 

5. Conclusions 

We have studied superconductivity under uniaxial compression in three /3-type salts, 
/3-(BEDT-TTF) 2 I 3 , /3-(BDA-TTP) 2 SbF 6 , and /3-(BDA-TTP) 2 AsF 6 , theoretically using 
both the two-band and dimer Hubbard models with FLEX approximation. 

We have found that the behavior of T c in /3-(BEDT-TTF) 2 I 3 with a stronger dimer- 
ization is reproduced by the dimer model, while that in weekly dimerized /3-BDA-TTP 
salts is reproduced by the two-band model. The present study will serve to clarify 
whether the dimerization of molecules is a prerequisite for high T c in organic supercon- 
ductors. 

Although we apply the uniaxial compression monotonically, the shift of T c exhibits 
a nonmonotonic behavior against the pressure. This result is theoretically considered 
to be due to the competition between the spin frustration and the effect of self-energy 
induced by the spin fluctuation. 
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